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' We evaluate ©(ojcts) corrections in the MSSM Higgs boson sector, including an analysis of the renormalization in 

^ , the bottom/scalar bottom sector. For jj, <0 the genuine two-loop corrections to the mass of the lightest Higgs boson 

mass can amount up to 3 GeV. Different renormalization schemes are applied and numerically compared. Likewise the 
' residual dependence on the renormalization scale is investigated. This allows to determine the remaining theoretical 

^ 1 uncertainties from unknown higher-order corrections at O(o(,a^) for different regions of the MSSM parameter space. 
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1. INTRODUCTION 
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' A crucial prediction of the Minimal Supersymmctric Standard Model (MSSM) [1] is the existence of at least one 
' light Higgs boson. Direct searches at LEP have already ruled out a considerable fraction of the MSSM parameter 



space [2, 3], and the forthcoming high-energy experiments at the Tevatron, the LHC, and the International Linear 



i-S^ ' Collider (ILC) will either discover a light Higgs boson or rule out Supersymmetry (SUSY) as a viable theory for 
^p^' physics at the weak scale. Furthermore, if one or more Higgs bosons are discovered, bounds on their masses and 
^-1. couplings will be set at the LHC [4, 5]. Eventually the masses and couplings will be determined with high accuracy 
(— I [ at the ILC [6-8]. Thus, a precise knowledge of the dependence of masses and mixing angles in the MSSM Higgs 
^ ■ sector on the relevant supersymmetric parameters is of utmost importance to reliably compare the predictions of the 
MSSM with the (present and future) experimental results. Even nowadays a precise prediction of the lightest MSSM 
Higgs boson mass is important to set solid bounds on the MSSM parameter space [9, 10]. 



The status of the available results for the higher-order contributions to the neutral CT'-even MSSM Higgs boson 
masses has been summarized in Refs. [11-13]. In particular two-loop contributions to the leading one-loop corrections 
involving the top and bottom Yukawa couplings (with at = hHiAn), ai, = /i^/(47r), ht^i, being the superpotential top 
or bottom coupling) have been evaluated. Corrections from the bottom/sbottom sector can also give large effects, in 
particular for large values of tan/3, the ratio of the two vacuum expectation values, tan/3 = V2/V1, and large values 
of ji, the supersymmetric Higgs mass parameter. The leading strong corrections at 0{abas) have been derived in 
Ref. [14] (in the limit tan/3 00) and in Ref. [15] (for arbitrary tan/3). In the (s)bottom corrections the all-order 
resummation of the tan /3-enhanced terms, C(af,(as tan/3)"), is also performed [16, 17]. Furthermore a full effective 
potential two-loop calculation exists [18], however, no public code is available. 

The potentially large size of corrections from the b/b sector makes it desirable to investigate the corresponding 
two-loop corrections, including the applied renormalization. An inconvenient scheme can give rise to artificially 
large corrections, whereas a convenient scheme absorbs the dominant contributions into the one-loop result, and 
higher-order corrections remain small. The comparison of different schemes (without artificially enhanced correc- 
tions) and the renormalization scheme dependence give an indication of the possible size of missing higher-order 
terms of 0{abal). 
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2. THE MSSM HIGGS SECTOR AND RENORMALIZATION 



The Higgs sector of the MSSM [19] comprises two neutral CT'-even Higgs bosons, h and H {rriii < tuh), the CT'-odd 
A boson (throughout this paper we assume that CV is conserved) and two charged Higgs bosons, H^. 

In the Fcynman-diagrammatic (FD) approach, the higher-order corrected Higgs boson masses, Mh and AIh, are 
derived as the poles of the /i, iJ-propagator matrix, i.e. by solving the equation 



-•hH 



0. 



(1) 



The renormalized self-energies, S^, can be expanded according to the one-, two-, . . . loop-order contributions, 

= E(i)(p2) + s(2)(p2)_^... ^s = hh,hH,HH . (2) 

The leading two-loop corrections from the b/b sector are the C(as) corrections to the dominant one-loop contributions 
oiO{ab)- They are obtained for zero external momentum and neglecting the gauge couplings, see Ref. [15] for further 
details. This approach is analogous to the way the leading one- and two-loop contributions in the top/stop sector 
have been obtained (see e.g. Rcf. [20]). Details about the renormalization in the MSSM Higgs sector can be found 
in Ref. [15]. 
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Figure 1: Some generic two- loop diagrams for the Higgs-boson self-energies 
arise for the tadpole contributions. 
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Figure 2: Some generic one-loop diagrams with counterterm insertion for the Higgs-boson self-energies {(j> ~ h, H, A), i,j, k = 
1,2). Similar diagrams arise for the tadpole contributions. 

The genuine two- loop Feynman diagrams to be evaluated for the Higgs boson self-energies, ^hh,hH,HHi are shown 
in Fig. 1. The diagrams with subloop renormalization are depicted in Fig. 2. The diagrams and the corresponding 
amplitudes have been generated with the package FeynArts [21, 22]. The further evaluation has been done using the 
program TwoCalc [23]. 

We now turn to the renormalization in the scalar-quark sector, entering via the diagrams depicted in Fig. 2. Since 
the two-loop self-energies arc evaluated at O^at^bCts) it is sufficient to determine the countcrtcrms induced by the 
strong interaction only. The tree-level relations and further details about our notation can be found in Rcf. [15]. In 
the following rriq., i = 1,2 denotes the two mass eigenvalues of the stop {q = t) and sbottom {q = b) sector. The 
corresponding unitary diagonalization matrix is denoted as Uq-j, i,j = 1,2. 

2.1. Renormalization of the top and scalar-top sector 

The t/t sector contains four independent parameters: the top-quark mass nit, the stop masses m^^ and m^^, and 
either the squark mixing angle 9^ or, equivalently, the trilinear coupling At- Accordingly, the renormalization of 



PSN0102 



this sector is performed by introducing four countcrtcrms that are determined by four independent renormalization 
conditions. The foUowing renormahzation conditions arc imposed (sec also Rcf. [24]). 

(i) On-sheU renormahzation of the top-quark mass 

(ii) On-sheh renormahzation of the stop masses 



(iii) The countcrterm for the mixing angle, 



ReEr (mH+ReSr (m? ) / 2(m? - m? ) 



Having already specified 5di, the At counterterm cannot be defined independently but follows from the relation 
sin 26*4 = [2mt{At - /iCot/3)] /[m? - m? ] . 



2.2. Renormalization of the bottom and scalar-bottom sector 

Because of S'?7(2)-invariance the soft-breaking parameters for the left-handed stops and sbottoms are identical, 
and thus the t and h masses are not independent but connected via the relation 

cos^ dirn^^ + sin^ ^i^^f, — cos^ + sin^ ^t*^^, + ^1 ~ ~ ^'^w cos(2/3) . (3) 

Since the t masses have already been renormalized on-shell, only one of the b mass countcrtcrms can be determined 
independently. In the following, the 62 mass is chosen as the pole mass, yielding Sm? ~ Rel]g^^(m? ), whereas the 
counterterm for to^^ is determined as a combination of other countcrtcrms. according to 

6mi = — \ — (cos^ dfSm- + sin^ 0fSmf - s'm^ OrSml - sin 26lr(m~ - nif )S6; 

61 cos Oij\ ^ tl I f2 02 t\ tl t2' ^ 

+ sm26f^{m^^ - mf^JSO^^ - 2m.t5mt + 2mbSmb^ . (4) 

Consequently, the numerical value of m^;^ does not correspond to the pole mass. The pole mass can be obtained from 
TOg^ via a finite shift of 0{as) (see e.g. Ref. [25]). 

There are three more parameters with counterterms to be determined: the fe-quark mass m^, the mixing angle 9^^, 
and the trilinear coupling Ai,. They are connected via 

sin26l^ = [2mbiAb ~ fitan(3)] l{m\^ - m?J , (5) 

which reads in terms of counterterms 

2 cos 26^7 06^7 = sm 26^7 1 ^ ^ sm26'r 7^ 5-^. (6) 

mb m~ — ?n~ mf — mf 

61 62 Oi 62 

Only two of the three counterterms. Smb, 50ij^, 5Ab can be treated as independent, which offers a variety of choices. 
In the following, four different renormalization schemes, collected in Tab. I, will be investigated. Details about the 
four schemes can be found in Rcf. [15]. 



Table I: Summary of the four renormalization schemes for the bottom quark/squark sector. Blank entries indicate dependent 
quantities. For more details see Ref. [15]. 



scheme 


62 




Ab 


e-b 


analogous to t/t sector ("mt OS") 


on-shell 


on-shell 




on-shell 


DR bottom-quark mass ( "mt DR" ) 


on-shell 


DR 


DR 




DR mixing angle and At ("Ai,, 6^ DR") 


on-shell 




DR 


DR 


on-shell mixing angle and Ab ("^6, OS") 


on-shell 




on-shell 


on-shell 
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2.3. Resummation in the 6/6 sector 



The relation between the bottom-quark mass and the Yukawa couphng hi,, which in lowest order reads mf, = 
/ifjUi/v^, receives radiative corrections proportional to hbV2 = hbta.n/3vi. Thus, large tan /?-cnhanccd contribu- 
tions can occur, which need to be properly taken into account. As shown in Refs. [16, 17] the leading terms of 
0{ab{as tan/3)") can be resummed by using an appropriate effective bottom Yukawa coupling. 

Accordingly, an effective bottom-quark mass, m|^^''^^^'^, is obtained by extracting its UV-finite tan /3-enhanced 
contribution Amt, (which enters through Shg) and writing it as 1/(1 + Am;,) into the denominator. In this way the 
leading powers of (a^tan/?)" are correctly resummed [16, 17]. This yields 

^DR,MSSM(^DR) ^ [^Ms (M^)5shift + l™,(sfin(^^2) ^ ^^^^(m,^)) + E^,"(m^)]/[1 + Am,] , (7) 
where = E^g + Am, denotes the non-enhanced remainder of the scalar 5-quark self-energy at 0(0;^), and b^^^^^ 



is given in by ft^'hift 
given at 0{as) by 



I - In 



3 Mi 



The tan /3-enhanccd scalar part of the 6-quark self-energy, Amt,, is 



m? TO? l0g(TO? /to? ) + TO? TO?log(TO? /to?) -I- TO?TO? log(TO?/TO? ) 



Am,^^2/i37T)aA^nPf,ms^^ , I \% / ^ ' . (8) 

TO- — mi){mi — m~ )(to- —to- I 

We incorporate the effective bottom-quark mass of (7) into our one-loop results for the rcnormalized Higgs boson 
self-energies. In this way the leading effects of 0{abas) are absorbed into the one-loop result. We refer to the genuine 
two- loop contributions, which go beyond this improved one-loop result, as "subleading 0{abas) corrections" in the 
following. 



3. NUMERICAL RESULTS 

For our numerical analysis we use the following parameters (if not indicated differently): rrit = 174.3 GcV, to^^ = 
2.94 GeV, Ma = 700 GeV, ^l = -1000, tan/3 = 50, MgusY = 1000 GeV, Af = 2Msusy, M2 = 2Mi = 100 GcV, 
rrig — 1000 GeV, fjP^ = mt. Large values of tan/3 and \^\ are chosen in order to illustrate possibly large effects in 
the b/h sector. The inclusion of all known corrections and the current experimental central top quark mass value of 
TOt = 178.0 GeV in our analysis would yield an increase in Mh of C(8 GeV) [11]. Therefore the mass values given in 
our numerical analysis should not be viewed as predictions of Mh] they are rather illustrations of the ag-corrections 
to the bottom Yukawa contributions at the two-loop level. 

In Rcf. [15] it has been shown that the "to, OS" scheme should be discarded. The reason for the problematic 
behavior of this scheme is easy to understand. The renormalization condition in the "jTifc OS" scheme is a condition 
on the sbottom mixing angle 9^ and thus on the combination (Af, — /ztan/3). In parameter regions where ^tan/3 is 
much larger than Ab, the counterterm 5Ab receives a very large finite shift when calculated from the counterterm 
59j^. This problem is avoided in the other renormalization schemes introduced in Tab. I, where the renormalization 
condition is applied directly to Ab, rather than deriving 5Ab from the renormalization of the mixing angle. 

In Fig. 3 it can be seen that the behavior of the subleading corrections in the three remaining renormalization 
schemes strongly depends on the choice of rrig. For rrig ^ 1000 GeV all schemes lead to an increase of Mh from the 
subleading 0{abas) corrections. For rrig > 1500 GeV, on the other hand, all schemes lead to a decrease, where the 
size of the individual corrections also strongly varies with mg. Accordingly, the relative size of the corrections in 
the different schemes also varies with rrig. Corrections up to about 3 GeV are possible. The differences between the 
three schemes are of 0(2 GcV) for large rrig. It should be noted that the effects of the higher-order corrections to 
Mh do not decouple with large rrig. The corrections at 0{atas) [20] as well as ©(a^as) grow logarithmically in the 
renormalization schemes that we have adopted. 

The jjP^ variation is shown in Fig. 4. The leading contribution (the 0{atas) result including resummation) is 
shown as the dark shaded (black) band. The results including the subleading corrections in the "mf, DR" scheme are 
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Figure 3: Mh as a function of rrig. 



shown as a light shaded (red) band. It can be seen that the variation with ij.^^ is strongly reduced by the inclusion 
of the subleading contributions. The variation with within the "m^ DR" scheme is tiny for rrig < 500 GeV, and 
reaches ±2 GeV for large rrig values. Thus, the fi^^ variation causes a similar shift in Mh as the comparison between 
the three renormalization schemes discussed above. 

The comparison of the results in the different schemes that we have analyzed and the investigation of the renormal- 
ization scale dependence give an indication of the possible size of missing higher-order corrections in the b/b sector. 
For ^ > (see Ref. [15]) the higher-order corrections from the b/b sector (beyond 0{abas)) appear to be sufficiently 
well under control. For /x < 0, on the other hand, sizable higher-order corrections from the b/b sector are possible. 
The size of the individual corrections and also the difference between the analyzed schemes varies significantly with 
the relevant parameters, /x, tan/3, uig and Ma- Wc estimate the uncertainty from missing higher-order corrections 
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Figure 4: dependence of Mh as a function of rrig. The black area corresponds to the 0{atas) result including resummation, 
i.e. the result without the subleading two-loop 0{abas) terms. The parameters are the same as in Fig. 3. 
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in the b/b sector to be about 2 GeV in this region of parameter space. 

The resuhs obtained will be implemented into the Fortran code FeynHiggs [11, 20, 26]. 
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